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Abstract 



We search for the Ml radiative transition 0(3686) — > jr] c (2S) by reconstructing the exclu- 
sive r] c (2S) ->• ^J^ ± 7r =F 7r + vr~ decay using 1.06 x 10 s 0(3686) events collected with the BESIII 
detector. The signal is observed with a statistical significance of greater than 4 standard de- 
viations. The measured mass of the r] c (2S) is 3646.9 ± 1.6(stat) ± 3.6(syst) MeV/c 2 , and the 
width is 9.9 ± 4.8(stat) ± 2.9(syst) MeV/c 2 . The product branching fraction is measured to 
be £(V>(3686) 7^(25")) x B(r/ C (25) K^K^w+n-) = (7.03 ± 2.10(stat) ± 0.70(syst)) x 
10 -6 . This measurement complements a previous BESIII measurement of ■0(3686) ^r] c {2S) with 
n c (2S) K^K^ and K+K'tt . 

no PACS numbers: 13.20.Gd, 13.25.Gv, 14.40.Pq 
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in I. INTRODUCTION 



112 Compared to other charmonium states with masses below the open charm threshold, 

n3 the properties of the rj c (2S) are not well-established. The determination of the rj c (2S) 

H4 mass, in particular, provides useful information about the spin-spin part of the charmonium 

us potential. The r) c (2S) was first observed at i?-factories |l|-|4j and, to date, the only two 

lie measured branching fractions are for decays to KKn and K + K~7i + n~n° [5|. While the 

ii7 absolute branching fractions currently have poor precision, BaBar used the two-photon 

us fusion process to measure the ratio of B(rj c (2S) — > K + K~ ir + 7r~7r°) to B(rj c (2S) — > K^I^n 11 ) 

n9 to be 2.2 ± 0.5(stat) ± 0.5(syst) Jdj. The production of the rj c (2S) is also expected from 

120 magnetic dipole (Ml) transitions |7|] of the ^(3686), and -^(3686) — > ^r] c {2S) with rj c (2S) — > 

121 KKtt has previously been observed by BESIII This analysis complements the previous 

122 analysis by focusing on the same radiative decay, ^(3686) — > jtj^S), but with r) c (2S) — > 

123 K^K+TT^TT+ir- . 

124 In our study, ^(3686) mesons are produced by the annihilation of electron-positron pairs 

125 at a center-of-mass energy of 3686 MeV. The production of the r] c (2S) through a radiative 

126 transition from the ^(3686) requires a charmed-quark spin-flip and, thus, proceeds via a 

127 Ml transition. Some of the generated t] c (2S) mesons will decay into hadrons, and then 

128 ultimately into detectable particles, like pions, kaons, and photons. We study the decay 

129 exclusively by reconstructing the r] c (2S) from its hadronic decay products and analyze the 

130 T] C (2S) candidate mass for an evidence of ^(3686) — > 7?7 C (2S'). The experimental challenge 

131 of the measurement of this decay channel is to detect the 48 MeV radiative photons in an 

132 experimental environment with considerable backgrounds, therefore the success of this study 

133 depends on a careful and detailed analysis of all possible background sources. 



134 II. THE EXPERIMENT AND DATA SETS 



135 The data sample for this analysis consists of 1.06 x 10 8 events produced at the peak of 

136 the ^(3686) resonance [9|. Data were collected with an additional integrated luminosity of 

137 42 pb~ x at a center-of-mass energy of ^=3.65 GeV to determine non-resonant continuum 

138 background contributions. The data were accumulated with the BESIII detector operated 

139 at the BEPCII e + e" collider. 



140 



142 
143 



The BESIII detector, described in detail in Ref. [10], has an effective geometrical accep- 
141 tance of 93% of An. It contains a small cell helium-based main drift chamber (MDC) which 
provides momentum measurements of charged particles; a time-of-flight system (TOF) based 
on plastic scintillator which helps to identify charged particles; an electromagnetic calorime- 

144 ter (EMC) made of Csl (Tl) crystals which is used to measure the energies of photons 

145 and provide trigger signals; and a muon system (MUC) made of Resistive Plate Chambers 
we (RPC). The momentum resolution of the charged particles is 0.5% at 1 GeV/c in a 1 Tesla 

147 magnetic field. The energy loss (dE/dx) measurement provided by the MDC has a reso- 

148 lution better than 6% for electrons from Bhabha scattering. The photon energy resolution 

149 can reach 2.5% (5%) at 1 GeV in the barrel (endcaps) of the EMC. And the time resolution 

150 of the TOF is 80 ps in the barrel and 110 ps in the endcaps. 

151 Monte Carlo (MC) simulated events are used to determine the detection efficiency, op- 

152 timize the selection criteria, and study the possible backgrounds. The simulation of the 

153 BESIII detector is based on geant4 [Uj , in which the interactions of the particles with the 

154 detector material are simulated. The ^(3686) resonance is produced with kkmc [l2[, which 
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is the event generator based on precise predictions of the Electroweak Standard Model for 
the process e + e~ — > // + 717, where / = e, /i, r, d, u, s, c, b, and n is an integer number. 



The subsequent decays are generated with EvtGen [131 ] . The study of the background is 
based on a sample of 10 8 ^(3686) inclusive decays, generated with known branching fractions 
taken from the Particle Data Group (PDG) [Bj], or with lundcharm for the unmeasured 
decays. 



lei III. EVENT SELECTION 

162 The decays of -^(3686) — > •yr] c (2S) with r) c (2S) — > if <jif ± 7r =F 7r + 7r _ are selected for this 

163 analysis. A charged track should have good quality in the track fitting and be within the 

164 angle coverage of the MDC, | cos 6*| < 0.93. A good charged track (excluding those from 

165 Kg decays) is required to pass within 1 cm of the e + e~ annihilation interaction point (IP) 
we in the transverse direction to the beam line and within 10 cm of the IP along the beam 
167 axis. Charged-particle identification (PID) is based on combining the dE/dx and TOF 
lee information to the variable yl m (i) = ( dE l dx ™°*^- dE l dx ^^ A ^ + / TOF measured -TOF expectcd ^ 

169 The values Xpid(0 an d the corresponding confidence levels ProbpiD(i) are calculated for 

no each charged track for each particle hypothesis % (pion, kaon, or proton). 

171 Photon candidates are required to have energy greater than 25 MeV in the EMC both 

172 for the barrel region (| cos^l < 0.8) and the endcap region (0.86 < | cos6*| < 0.92). In order 

173 to improve the reconstruction efficiency and the energy resolution, the energy deposited in 

174 the nearby TOF counter is included. EMC timing requirements are used to suppress noise 

175 and remove energy deposits unrelated to the event. Candidate events must have exactly six 

176 charged tracks with net charge zero and at least one good photon. 

177 K s candidates are reconstructed from secondary vertex fits to all the oppositely charged- 

178 track pairs in an event (assuming the tracks to be 7r =t ) . The combination with the best fit 

179 quality is kept for further analysis, where the K s candidate must have an invariant mass 

180 within 10 MeV/c 2 of the K s nominal mass and the secondary vertex is well separated from 

181 the interaction point. At least one good K s is reconstructed, and the related information is 
used as input for the subsequent kinematic fit. 

After tagging the ir + 7r~ pair from the K s , the other charged particles should be three 

184 pions and one kaon. To decide the species of those particles, we make four different particle 

185 combination assumptions: K + n~n + Ti~ , 7t + K~7t + tt~ , tt + tt~ K + tt~ , and tt + 7t~tt + K ~. For 
lee the different assumptions, four-momentum conservation constraints (4C) are required to 

187 be satisfied for each event candidate. For each event, the Ml-photon is selected with the 

188 minimum chi-square of the 4C kinematic fit (xlc) by looping over all the good photons. Then 

189 the xlc an d th e chi-squares of the particle-identification for kaon (Xk) an d pions (% 2 ) are 

190 added together as the total chi-square (Xtotai) f° r even t selection. The types of particles are 

191 determined by choosing the smallest total chi-square. Events with Xtotai < 60 are accepted 

192 as the / yK s K ± 7r T 7T + 7T~ candidates. 

193 To suppress the ^(3686) — > n + n~ J/ip, J/ip — >■ r yK s K ± n T decay, events are rejected 

194 if the recoil mass of any ti + tx~ pair is within 15 MeV/c 2 of the J/ip nominal mass. The 

195 ^(3686) — > rjj/ip, 7] — > 77r + 7r~ events are rejected if the mass of K s K ± n T is greater than 
we 3.05 GeV/c 2 . In order to suppress ^(3686) — > r]'K s K ± ir T , 77' — > 77r + /T~ decays, events are 
197 removed if the mass of any 77r + /T~ combination is within 20 MeV/c 2 of the nominal 7/ mass. 
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198 IV. DATA ANALYSIS 



199 The results of an analysis of the inclusive MC data sample showed that the primary 

200 source of background is ^(3686) — > KgK ± 7i T 7i + 7i^ . There are two mechanisms for this de- 

201 cay to produce background: a fake photon, or a photon from final-state radiation (FSR) is 

202 incorporated into the final state. Other backgrounds include ^(3686) — > 7r° KgK ±/ K T 7i + 7i^ 

203 with a missing photon and initial state radiation (ISR). The phase space process ^(3686) — > 

204 r yKgK ± 7i T 7Y + n~ has the same final states as our signal, so it should be considered as an ir- 

205 reducible background. As discussed in a later section, the size of this irreducible background 

206 is estimated using a region of KgK ± ir T 7i + 7i^ mass away from the rj c (2S) mass. 

207 In the ^(3686) — > KgK ± -K T ir + TT~ background with a fake photon, a peak could be pro- 
2ob duced in the K® K ± tt t it + it~ mass spectrum close to the expected r] c (2S) mass with a sharp 

209 cutoff due to the 25 MeV photon energy threshold. Considering that the fake photon does 

210 not contribute useful information to the kinematic fit, we set the photon energy free in the 
2n kinematic fit to avoid the mass distortion caused by the 25 MeV photon energy threshold. 

212 We call this the 3C kinematic fit and produce the mass spectrum based on it. MC studies 

213 demonstrate that with the 3C kinematic fit, the energy of the fake photon tends to zero, 

214 which is helpful in separating the signal from the fake photon background, as shown in 

Fig mM- 
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FIG. 1: Invariant mass spectrum of K^K^-ir^^ir" for the background ^(3686) — > KgK 7r =F 7r + 7r~ 
with a fake photon (left panel) and the signal ^(3686) — > jr] c (2S), rj c (2S) — > K^K^^^tt' (right 
panel). The points with error bars are 3C kinematic fit results, and the solid lines are 4C kinematic 
fit results. 

In the other ^(3686) — > KgK ± 7r T 7r + 7r~ background, a photon from final state radi- 
ation (7fsr) could contaminate our signal. The M^ K3n with the FSR process has a 

long tail from 3.58 GeV/c 2 to 3.68 GeV/c 2 in our r) c (2S) signal region. We have to es- 
timate the contribution of this FSR process, because it contributes to the background in 
our signal region and cannot be reduced for the same final states as the signal. FSR is 



221 simulated in our MC generated data with PHOTOS [15] . and the FSR contribution is 

222 scaled by the ratio of FSR fractions in data and MC generated data for a control sample 

223 of ^(3686) — > / -fK + n~K + K~ and ^(3686) — > 77r + 7r~7r + 7r~ [l6[. The background contribu- 

224 tions from ^(3686) — > K^K ± h t 'k + 'k~ with fake photons and 7fsr are estimated with MC 

225 distributions normalized according to branching ratios we measured. 
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226 The channel ^(3686) 



227 the photons from the 



228 7T 



7T _ can contaminate our signal when one of 
7T° is not detected. MC generated events of the -^(3686) — > 
KgK ± -ir T ii + 7i^ process, based on the phase space model, and which satisfy the selec- 
tion criteria for the ^(3686) — > 'yKgK ± 7r T 7r + 7T~ signal, are taken to study this background 
and estimate its response. To prove the correctness of the MC simulation, the ^(3686) — > 
ir°KgK ± 7T T 7T + iT~ control sample, which is selected from the colliding data, times the effi- 
ciency to reconstruct ^(3686) — > 7r KgK ± ii T 7i + ii~ events as ^(3686) — > r yKgK ± TT T 7i + ir~ is 
shown in Fig. [2] and compared with the same distribution obtained from the corresponding 
^(3686) — > -K°KgK ± ir T 7i + 7i^ MC simulation. The consistency of the two distributions is 



checked by the Kolmogorov-Smirnov test [17J, and a good agreement is verified (the consis- 
tency probability reaches 0.28). 



i i i i i i i i i i i i i i i i i i i i i i i i i i i i i 




K£K3tc 
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FIG. 2: The invariant mass distribution of K^KSir for the background from t/>(3686) — > 
. The black circles with error bars show the background shape obtained from 
the collider data. The red triangles with error bars represent the M K o K3lT distribution from a 
corresponding MC sample. 

237 The background from the continuum (including ISR) is estimated with collider data 

238 taken at a center of mass energy of 3.65 GeV. The events must pass the signal selection 

239 requirements and are then normalized according to differences in integrated luminosity and 

240 cross section. Particle momenta and energies are scaled to account for the beam-energy 

241 difference. The resultant number and the KgK3ir invariant mass shape considering these 

242 scale factors (/continuum = 3.6) are used in the final fit. 

243 The background from phase space has the same final states as the signal. To select a 

244 clean phase space sample, the M K o K3w region [3.20,3.30] GeV/c 2 is chosen. This choice is 

245 made because there is a long r\ c tail in the area M K o K37r < 3.0 GeV/c 2 which originates from 

246 the decay channel -^(3686) — > jric. There are three obvious peaks in the area M K o K3n > 
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247 3.3 GeV/c 2 which are from the decay channel ^(3686) — > 'yXcJ, {J = 0,1, and 2). The 

248 branching fraction of the phase space process is calculated to be 1.73 x 10" 4 . The KgK3ir 

249 invariant mass spectrum of MC phase space events is used in the final fit, while the number 

250 of events is left floating. The number of phase space events obtained by fitting the mass 

251 spectrum is consistent with that estimated by the branching fraction we calculated. 

252 In the KgK37i mass spectrum fitting, the fitting range is from 3.30 GeV/c 2 to 3.70 GeV/c 2 

253 so that the contributions of backgrounds and Xcj{J = 0,1, and 2) can be taken into account. 

254 The final mass spectrum and the fitting results are shown in Fig. [3j The fitting function 

255 consists of the following components: r] c (2S), Xcj(J = 0,1, and 2) signals and ^(3686) — > 
KgK ± ir T 7T + 7r~ , ^(3686) — > 7r°KgK ± iT T 7r + 'K~' , ISR, and phase space backgrounds. The line 




FIG. 3: The results of fitting the mass spectrum for Xcj and rj c (2S). The black dots are the col- 
lider data, the blue long-dashed line shows the Xcj and rj c (2S) signal shapes, the cyan dotted line 
represents the phase space contribution, the violet dash-dotted line shows the continuum data con- 
tribution, the green dash-double-dotted line shows the contribution of ^(3686) — > KgK ± 7r T n + ir~ , 
and the red dashed line is the contribution of ^(3686) — > ir° "K^K^ir^ : ?r + 'tt~ . 

256 

257 shapes for Xcj are obtained from MC simulations. These can describe the Xcj spectrum 



well in the collider data after applying the kinematic fit correction [18[. The line shape for 
i] c {2S) produced by such a Ml transition is given by: 

(E^ x BW(m) x damping (Ej)) <g) Gauss(0,a), (1) 



260 where BW(m) is the Breit-Wigner function, m is the invariant mass of KgK37r, E 1 = 

9 9 
7TL — TTl 

261 2m^ 8 ( 36 86) ^ s ^ ne ener gy °f the transition photon in the rest frame of ^(3686), damping(Ej) 

262 is the function to damp the diverging tail raised by E^ and Gauss(0, a) is the Gaussian 

263 function describing the detector resolution. The detector resolution is determined by the 



S 



MC study, and the difference of data and MC has been taken into account which introduces 
negligible uncertainties in branching fraction, mass and width measurements comparing with 
other factors. The form of the damping function is somewhat arbitrary, and one suitable 



function used by KEDR [19J for a similar process is 



E 2 

damping^) = EEq + { ° Eq)2 , (2) 



268 where E = \ — V ^ 2S) is the peaking energy of the transition photon. Another damping 



269 function used by CLEO [20j is inspired by the overlap of wave functions 

damping(E 7 ) = exp(—E 2 /8(3 2 ), (3) 

270 with = (65.0 ± 2.5) MeV from CLEO's fit. In our analysis, the KEDR function (Eq. |2J) is 

271 used in the fitting to give the final results, and the CLEO one (Eq. [3]) is used to estimate 

272 the possible uncertainty caused by the form of damping functions. 

273 The result for the yield of rj c (2S) events is 57 ± 17 with a significance of 4.2a. The 

274 significance is calculated from log-likelihood differences between fits with and without the 

275 r) c (2S) component. The robustness of this result was tested by considering different damping 

276 factor forms, FSR fractions, and background assumptions. In all the cases, the statistical 

277 significance is found to be larger than 4a. The resulting mass and width from the fit are 
27B 3646.9 ± 1.6 MeV/c 2 and 9.9 ± 4.8 MeV/c 2 (statistical errors only), respectively. We find 

279 the product branching fraction i3(^(3686) ->■ JVc(2S)) x B(r] c (2S) -> K^K^n+ir) = 

280 (7.03 ± 2.10) x 10~ 6 with the efficiency of 11.1% for the signal selection. 



281 V. ESTIMATION OF SYSTEMATIC UNCERTAINTIES 

282 The systematic uncertainties in the rj c (2S) mass and width measurements are estimated 

283 by the uncertainties in the damping factor, scale factor and the number of ^(3686) — > 

284 7T°KgK ± 'K T 7r + 'iT~ events. The results are summarized in Table [U and described in more 

285 detail in the following. 



TABLE I: Uncertainties in the mass and width of r] c (2S). 



Source 


mass uncertainty width uncertainty 


Damping factor 


< 0.1% 


28% 


Scale factor 


negligible 


5% 


No. of n°K%K ^n+w- 


< 0.1% 


5% 


Total 


< 0.1% 


29% 



286 We change the damping factor to the CLEO form, then compare the results with that 

287 obtained with the KEDR form, and the difference is taken as the uncertainty originating from 

288 the damping factor. The background shape of ^(3686) — > KgK ± n T n + n~ could influence the 

289 fitting results, so we change the FSR scale factor of 1.46 by la to 1.412 and 1.504, and the 

290 difference in the results is taken as the uncertainty coming from scale factor. In the fitting 
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291 of the mass spectrum, the number of events for if) (3686) — » TT°Kg K ± tt t 7t + ti^ is fixed. We 

292 change the number of events by la, and take the difference in the results as the uncertainty 

293 originating from the number of background events from ^(3686) — > 7r°K s K ± 7T T ii + 7i^ events. 

294 The systematic errors in the measurement of the branching fraction are summarized in 

295 Table [11] and explained below. 

TABLE II: Summary of systematic uncertainties in the measurement of 
£(V>(3686) -». 7Tfc(2S), r) c (2S) -»• K^K^Vtt') . 



Sources 


Systematic uncertainties 


MDC tracking 


4% 


Photon reconstruction 


1% 


K s reconstruction 


4% 


Kinematic fitting and PID 


2% 


Total number of ^(3686) 


0.8% 


Damping factor 


2% 


Scale factor 


5% 


No. of ^(3686) -> ^K^K^tt+tt- 


2% 


rj c {2S) width 


3% 


Intermediate states 


5% 


Total 


10% 



296 The tracking efficiencies for K ± and 7r ± as functions of transverse momentum have been 

297 studied with the process J /if) — > K s K ± n T , K s — > 7r + 7r~ and ^(3686) — > n + 7i~J/ip, respec- 
29B tively. The efficiency difference between data and MC is 1% for each K^ track or ir^ track 



299 



2ll . |22| . So the uncertainty of the tracking efficiency is 4% for four charged tracks. The 

300 uncertainty of the two pions from Kg is not included here, because it is included in the K s 

301 uncertainty. 

302 The uncertainty due to photon reconstruction is 1% per photon j23[. This is determined 

303 from studies of photon detection efficiencies in the process J/if> — > p°7r°, p° — > 7r + 7r~ and 

304 7T° — > 77. 

305 Three parts contribute to the efficiency for Kg reconstruction: the geometric acceptance, 

306 tracking efficiency and the efficiency of K s selection. The first part was estimated using 

307 an MC sample, and the other two were studied by the process J /if) — > K*K° + c.c. The 
3ob difference between data and MC is estimated to be 4%. 

309 To estimate the uncertainty of kinematic fitting, we first correct the track helix param- 

310 eters (<f>o, k, tg\) to reduce the difference on xlc fr° m kinematic fitting between data and 

311 MC, where 0o is the azimuthal angle specifies the pivot with respect to the helix center, k is 

312 the reciprocal of the transverse momentum and tgX is the slope of the track. The correction 

313 factors are obtained from J/ip -> 0/ o (98O), <p K + K~ and / (980) tt+tt". The MC 

314 samples after correction are used to estimate the efficiency and fit the invariant mass spec- 

315 trum. Fig. 0] (left) shows the xIc+pid distribution with and without the correction in MC 

316 and in data. The distribution of xIc+pid with correction is closer to the data than without 

317 correction. However, the agreement is not perfect, and we take the systematic uncertainty 



318 to be the difference of the efficiency between MC before and after correction [18J|. The com- 
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319 parison is shown in Fig. [4] (right). The systematic uncertainty from kinematic fitting is 2% 

320 with xIc+pid < 60. 




FIG. 4: [left panel] The xIc+pid distribution with and without the correction in MC and in data. 
The black dots show the distribution of xIc+pid m the data, the orange (green) histogram repre- 
sents the distributions of xIc+pid without (with) correction in MC. [right panel] Efficiency results 
with and without correction at different X^c+pid cu ts- 

321 We also change the form of the damping factor, the value of the FSR scale factor and 

322 the number of events for -^(3686) — > -K KgK ± Tc T ii + 7i^ to estimate the uncertainties in the 

323 branching fraction, which is the same as the method to estimate the uncertainties of r] c (2S) 

324 mass and width. The total number of ^(3686) events is estimated by the inclusive hadronic 

325 events, and the uncertainty is 0.8% jsf. 

326 To estimate the uncertainty due to the rj c (2S) width, we change the rj c (2S) width of 

327 9.9 MeV/c 2 by la to 5.1 MeV/c 2 and 14.7 MeV/c 2 in the MC simulation. Comparing 
32B the efficiencies with 11.1%, which is used in calculating the branching fraction, we find a 

329 difference of 3%. 

330 For the uncertainty from intermediate states, we generate MC samples including these 

331 states (if* (892), p) and compare the corresponding efficiencies. We take the 5% difference 

332 as the uncertainty. 

333 We assume that all the sources of systematic uncertainties are independent and the overall 

334 systematic uncertainties are obtained by adding all single ones in quadrature. 



335 VI. CONCLUSION 

336 We observe the decay mode r] c (2S) — > KgK ± 7i T n + 'K~ and establish the Ml transition 

337 of ^'(3686) — > •yrj c {2S) using this decay mode. The mass of the r] c (2S) is measured to be 
33B 3646.9 ± 1.6(stat) ± 3.6(syst) MeV/c 2 , and the width is 9.2 ± 4.8(stat) ± 2.9(syst) MeV. 

339 Comparing with BESIII previous measurements [si], the width is consistent with each other 

340 within 1 standard deviation and the mass is about 2 standard deviation. The product 

341 branching fraction is measured to be £>(-?/>(3686) — > jr] c (2S)) x B(rj c (2S) — > KgK ± 7i T 7i + n^) 

342 = (7.03 ± 2.10(stat) ± 0.70(syst)) x 10~ 6 . The statistical significance is greater than 4 

343 standard deviation. 
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344 To compare with the BABAR results 

B{r) c {2S) -)• K+K-7T+7T-7T ) 



B(r] c (2S) -> K° S K^) 



2.2 ±0.5 ±0.5, (4) 



34 5 we take the value of (4.31 ± 0.75) x 10" 6 as measured by BESIII for #0(3686) -> ir) c (2S)) x 

346 B{r} c (2S) -»■ iT^ ± 7T T ) 0, and assuming that 

B( Vc (2S) -» if+if-TT+Tr-Tr ) = 

i3(r/ c (25) -)• UrgA-iTrFTr+TT-) ' ' 1 J 

347 where the value 1.52 is calculated in Xcj decays, which has the same isospin, we obtain 

B( Vc (2S) -> _ g(„ c (2g) -> g^Wr) _ 

B(ifo(25) ^±ttF) " L52 BM25) ^ ^±ttT) " ^±0.5o±0.33. (6) 

34B These two results are consistent with each other after considering the statistical and sys- 
349 tematic uncertainties. 
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